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t:^- ! Abstract 

(N; 

r-| ' In a generalized Airy matrix model, a power p replaces the cubic term of 

'V ■ the Airy model introduced by Kontsevich. The parameter p corresponds to 

r^. Witten's spin index in the theory of intersection numbers of moduli space of 

^ ! curves. A continuation in p down to p = —2 yields a well studied unitary 

I matrix model. The application of duality and replica to the p-th Airy model 

CN ■ provides, through this equivalence, a generating function for both the weak 

Q ■ and the strong coupling expansions of the unitary model. We thereby recover 

cn . and extend further the results for these expansions. 

g '. 1 Introduction 

o: 

^ ■ The unitary matrix model with an external source is defined by the integral 

1> ■ over the unitary group U{N) 

'O ; 

"§: Z = |rff/e*^(^^'+^'^) (1.1) 

where U is a N x N unitary matrix, WU = 1. The external source matrix 
C is a given fixed complex matrix. Gross and Witten [1] have studied the 
case in which C is a multiple of the identity, C = N/g^, and showed that. 
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in the large N limit, there is a third order transition at some critical Qc- 
Brezin and Gross [2] generalized the model to include an arbitrary source 
matrix C and again found that there is a phase transition governed by the 
parameter tr(C^C)^^/^. This model has a weak coupling expansion , when 
the eigenvalues Aj of C^C are large, which involves powers and products of 
parameters 

t„. = ET;^ = Tr(C7tC)-(-V^). (1.2) 

i \ 

It has also a strong coupling expansion in powers of products of the param- 
eters Tr(C^C)™. The strong coupling expansion is straightforward and the 
coefficients involving powers of the A™ have been tabulated long ago, up to 
m = 5 [3]. The weak coupling expansion can be found by different methods ; 
Gross and Newman have used Virasoro equations [4] and we reproduce their 
analysis in section 4. We will follow here a completely different approach: 
we shall based the analysis on another matrix model, which will be shown to 
be equivalent to (1.1), whose interest is to show that the coefficients of the 
expansion of log Z in powers of the t^'s are related to topological invariants, 
namely the intersection numbers of the moduli space of p-spin curves on a 
Riemann surface, in the singular limit p — )■ —2. The same model will also be 
used to recover and extend the results of the strong coupling expansion. 

Indeed we shall demonstrate a remarkable equivalence of the unitary 
model (1.1), with a matrix model given by an integral over N x N Hermitian 
matrices 

Zb= I dBe''^^^''^''^^+''^^ (1.3) 

where B is a NxN Hermitian matrix, A = CC\ provided one tunes properly 
the parameter k. This equivalence was already discussed by Mironov et al. [5, 
6] through Virasoro equations and determinant formulations. We will show 
that one can recover it from an explicit determinantal expression in terms of 
Bessel functions [7], and from a comparison between the expansions based 
on (1.3) with that of (1.1). However the equivalence requires some proviso, 
as the reader, puzzled by the singularities due to vanishing eigenvalues of B, 
may guess. We are in fact dealing here with matrices which are unitarily 
equivalent to complex diagonal matrices. Therefore if we wrote B = UXU\ 
where X is diagonal with complex entries, the integral over the unitary group 
would yield a Jacobian which is the modulus square of the Vandermonde 
determinant J = |A(X)p = Y[i<j l^i ~ ^jl"^ and the integration over the XjS 
are over specified contours in the complex plane. 

Although the equivalence holds only if the coefficient k goes to zero in 
the weak coupling limit and to — A^ in the strong coupling limit, we keep it 
arbitrary at this stage. The weak coupling expansion for k = gives the 
intersection numbers of the moduli space, and for k ^ 0, the expansion is 



related to the intersection numbers of the discretized moduh space [8]. The 
expression of (1.3) is the special case, p = — 2 and r = —1, of a (p, r) model 



Z(„ ,) = / dBexpi —trBP+^ + — ^tr5^+^ + trfiA) (1.4) 

J p + I r + 1 



1 „„,i k 

P 

which generalizes Kontsevich model [9] p = 2,k = 0. The model that we 
want to study involves an analytic continuation from positive integer values 
of p and r to negative values. 

The reason that we find the model (1.3) interesting, in addition of the 
topological invariants that it allows one to compute, is that it leads to com- 
pletely different ways of calculating explicit expressions for the expansions. In 
our previous work we have discussed a powerful duality between expectation 
values of characteristic polynomials [10, 11, 12, 13], namely the relation: 

~ N „ n 

/ dM T] det(A„-M)e-^t'^^'+t'-^^ = / dB T] det(a.-iB)e-'^'''''+''''^ 

Jnxn ^J-^ JnxN jJ^ 

(1.5) 

where A = diag(Ai, ..., X^) and A = diag(ai, ..., a„). Both sides are integrals 
over Hermitian matrices (properly normalized). The left-hand side is an 
A^ point function for n x n matrices in an external given matrix source A, 
whereas the r.h.s. is an n-point function for N x N matrices, in the matrix 
source A. The identity (1.5) generalizes the relation between Hermitian and 
Kontsevich- Penner models aX A = [14, 15]. The strategy that we use 
consists of the following successive steps: 

• one tunes the numbers Oj in the r.h.s. in order to generate, in an 
appropriate large n scaling limit, the {p,r) model (1.4). 

• in the Gaussian model of the left-hand side we know explicit exact 
expressions for the correlators [10, 16] 

f/(si, ..., si) =< tre^^^tre^^^ ■ ■ ■ tre^'^ > . (1.6) 

in which the expectation value is meant with the normalised weight 
proportional to e~'2^^^ +trMA_ 

• the replica method, which consists of repeating the eigenvalue A^, n^ 
times and to let n^ go to zero. This is used to relate the characteristic 
polynomials of the l.h.s. to the usual correlation functions as in 

lim4-^(det(A-M))^ = tr— i— (1.7) 

The resolvent is related by Fourier transform to tre^^ which enters into 
(1.6). We will explain later the meaning of taking an A^ = limit here 
in spite of the fact that in the right-hand-side of (1.5) A^ is fixed. 



In this article, we use this strategy and tune the external source A to 
generate the(p, r) model of (1.4). The limit p = —2 and r = —1, obtained 
through a nontrivial continuation from positive values of p and r, is of par- 
ticular interest. This rather involved construction is finally compared with 
the known expansions of the unitary model Zu of (1.1) [3, 4, 5] both in the 
strong and in the weak coupling phases. Our results provide series expan- 
sion for intersection numbers with / marked points and arbitrary genus. The 
consistency with known results confirms that the duality plus replica method 
may be succesfully extended to negative values of p and r in the {p, r) model 
(1.4). 

2 Dual model and its /-point functions 

The starting point is the duality relation (1.5) 

„ N „ n 

/ dM n det(A„-M)e-^*=^^'+*'-*^^ = / dB H det(a,-2fi)e-^*^^'+*'-^^ 

(2.1) 
where A = diag(Ai, ..., X^) and A = diag(ai, ..., a„) ; the proof is simple and 
can be found in [10]. We now tune the aj in the r.h.s. to produce the model 
(1.4) . We use the identity 

n ■ n oo -m 

n det(l - -B) = exp(- E E ZZZ^^^B-) (2.2) 

j=l "i i=l m=0 "'"j 

and fix now the sums I]?=i ^- L^t us call p the first one(m = 1) 

p = ^J:-^ (2.3) 

which may be absorbed by a shift of A. Taking next Z]j=i ^ = 1 the term 

— ^tri?^ in the exponential is cancelled. Assume p < r; we can cancel all 
the terms other than tiB^ and tiB'^ by chosing the aj so that I]j=i ^ = 0, 

for all m < r, except m = 2 and m = p. Note that this requires in general 
complex ttj. In order to eliminate the terms with m larger than r, we have 
to consider the following large n scaling limit : 

• assume each Uj is of order n^/^ so that J2j -^ = 0(1) and in fact we 
chose them so that this last sum is exactly equal to one. 



we take the Aq, near p in a range in which A^ — p ~ n 



-l/2+2/r 



• we take the matrix B of order n^^^ ^/''. Then one has {J2 1/aptri?'' of 
order one 

• in this scale Y, l/a^'tiB"^ ~ n'^i^-™-/''") goes to zero for tti > r and one 
can drop all the powers higher than q. 

• however the term Yl 1/a^tri?^ ~ ^2(i-p/r) gj-Q-y^^rg -^ji;}^ ^^ since p < r and 

we have to tune the aj so that Y l/'^j instead of being of order n^^^'^ 
is of ordern^-P/2-2(i-P/r-)_ 

In this scaling limit, with the aj chosen as indicated, the r.h.s. gives the 
model Z(^p^r) of (1-4) • 

The method consists of using now the left hand side of (2.1). There the 
duality left us with the Gaussian model 



■'dual 



rfMe-^t^^'+t"-^^ (2.4) 



with the matrix A fixed as just discussed above. The reason for using this 
dual model is that this we have for (2.4), explicit expressions for the /-point 
correlation functions 

t/(si, ..., si) =< tre^i^tre^^^ ■ ■ ■ tre^'^ > (2.5) 

where the average stand for 



<X >= -^ fdMX{M)e- 

Zdii.nJ J 



itrM^+trMA 



■'dual 

The result, derived in [10, 11], is 

in which the contour integrals circle around all the poles aj. 
Let us first consider the one point function U{s). 



(2.6) 
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U{s) 



(2.7) 

We use now the aj tuned as discussed above. Expanding the logarithm, we 
have 



Using the previous tuning for the aj in the large-n scahng hmit we obtain 

— i (^^ .T^({u+sY+^-u^+^)+J^((u+sY+^-u-+^) 



For /-point function f/(si, ..., S/), we obtain similarly. 



U{si,...,si) 



J-r dUi _L 



/ /ii 227r 



,+s,)'-+i-«fi) 



X det . (2.9) 

Ui-Uk + Si 



We have studied earlier the case in which k = and p > ; there the 
/-point correlation functions U{si, ...,Si) are the generating function of the 
intersection numbers of the moduli space of curves with p-spins [12, 13]. 

3 Weak coupling expansion 

We now study the model in which one continues from integer positive p and 
r down to p = —2 and r = — 1. We have in this case, 

U{s) = -(f ^^-^s + l+klog{u+s)-k\ogu ^3^^^ 

In the original unitary model (1.1), a simple model of one plaquette lattice 

QCD, the matrix C is inversely proportional to the gauge coupling constant. 

Therefore the eigenvalues Aq of C'^C are large in weak coupling, small in 

strong coupling. The variables in U{s) is a Laplace transform of A; therefore 

the weak coupling limit corresponds to large A, hence small s. We first rescale 

U — 7- su. 

^ X / N f du 1 ,u + 1,1, , , 

U(s) = i e'--^-+'H^—f, 3.2 

J 2171 u 

since we continue later in k down to negative values, the contour integral is 
really an integral over the discontinuity along the real -u-axis. We chose a 
countour in the -u-plane which goes parallel to the imaginary axis through 
the point u = —1/2 at which l/u{u + 1) is maximum in the real direction 
and change variable -u = l/2(— 1 + i/x^/s) : 



1 f dx 43:^ 1 _ ix^/M. , 



This integral contains a singular term proportional to 1/v^ which we discard, 
i.e. in order to obtain the terms of the expansion with positive powers of s, 
we consider instead 



For s small, we can expand the integrand in powers of s and compute the suc- 
cessive coefficients as Gaussian integrals. The result for the positive powers 
of s takes the form 



1 

21-218' 



-s^—±—{Ae - 1)(4A;2 - 9)(4A;' - 25)(4A;2 - 49) 



~^' 135^^^^' - 1)(4A;2 - 9)(4A;2 - 25)(4A;2 - 49)(4A;2 - 81) + ©(s^) 

(3.5) 

If we set k = in this result the {p, r) model of (1.4) reduces to the p-th 
Airy model, with here p = —2. Let us verify that the coefficients of s™ in 
the above series (3.5) do coincide with the intersection numbers for p = —2 
of the p-th generalized Kontsevich model [13]. Indeed for this model we had 
found 

^ (p-I)(p^3)(1.2rt ^.„,_g,^..., ,3,, 

where y = s ~^p. For p = —2, this does agree with the series (3.5) for k = 0. 
In the next section we shall compare these results with the weak coupling 
expansion derived both for the model defined by (1.3) and for the unitary 
model (1.1) and see that they agree with (3.5) in the limit k = 0. The alleged 
correspondence of these two models is in fact true in the limit k — )■ — A^, but in 
this calculation we are focusing on the terms of the weak coupling expansion 
which contain only one single trace. The flow of internal index has only one 
cycle and in the zero replica limit only those terms survive. 



4 Exact calculation, Virasoro constraints and 
comparison 

Instead of using duality and replica we return to the partition function Z^ 
defined in (1.3). Standard techniques allow us to express Z^ as a determinant 
of modified Bessel functions for arbitrary k. We shall compare the weak 
coupling expansion deduced from this explicit expression, with that obtained 
above from f/(s) using the dual model. Denoting the eigenvalues of B by Xj 
(j = 1, ...N), we can integrate over the unitary group in the matrix integral 
(1.3). Given the non commutation of the matrices B and A this requires the 
use of the HarishChandra (Itzykson-Zuber) formula : there is thus a division 
of the Jacobian measure |A(x)p by A(a;) . Once the integration over the 
unitary group is done we have 




N ^TT — T N 

^B= /n^^^xSn^^^'*^^"''- (4-1) 

in which the complex eigenvalues Xi circle around the origin in the complex 
plane. Changing variables Xi — )■ -^e*^% the integrals over these angles pro- 
duce modified Bessel functions Im{z), and Zb has an explicit expression as 
a a determinant, 



'B = TTTTTl+rTTTT^ex I — rrJfc+il^V ""i) I (4-2) 

iUK) 2 A (A) \x: 

In the weak coupling region of large Xj, we use the asymptotic expansion of 
the modified Bessel functions: 

/^(2v^) = ^^(1 - —^ + ^- il^!—^ 

vWX 4v^ 2!(4v/A)^ 

3!(4v^)3 '' ^ ' ' 

Retaining only the leading terms in (4.3) for large Xj, gives the term of order 

N 1 N 



zero 

This is the nothing but the genus zero contribution to the free energy F 
found in [2] for the unitary model. 



Keeping N finite, we can calculate a few terms in the expansion in inverse 
powers of A: 

(2A; + 2iV)2 - 1 ^ 1 



16 ^v^ 

((2A; + 2iV)^-l)((2A; + 2iV)^-9) J^ 1 , 
^ 512 ^{r[ v^' 

{{2k + 2Nf - l){{2k + 2Nf - 9) 
^ 3!46 

X -8E77Tr + ((2^ + 2iV)'-17)(E 



+ ■ ■ ■] (4.5) 

Note that A^ appears always in the combination {k + N). Therefore, setting 
k = —N , all the factors (2A; + 2iV) vanish and we find that the weak coupling 
expansion of the partition function Zb (1-3), up to this order , agrees with 
the weak coupling expansion obtained above in (3.5) with the dual model at 
A; = 0, if we look only at the single trace terms. Furthermore it agrees with 
the known weak coupling expansion of the unitary model [7] Zu defined by 
(1.1). This is as a manifestation of the equivalence of Zu in (1.1) and Zb in 
(1.3) at A; = -A^. 

This same equivalence at A; = — A^ in (4.5) can be established through 
the weak coupling expansion based on Virasoro equations. The Virasoro 
equations for the partition function Zu in (1.1) have been derived by Gross 
and Newman [4]. We quote here their results; the unitary model satisfies the 
obvious equation 



d^Zr 



T- = 6acZu (4.6) 

but Zu is only function of the eigenvalues Aj of CC''' ; one thus obtains [2] 



d Zu y^ I dZu dZu ^ 1 .y v^ dZu . , . 

With Zu = ZqY, the Virasoro constraints [4] take the form 

L„r = -dnY n > (4.8) 

with the Ln given by 

°° 1 ~ 1 

n = : ^(A; + -)tAY + — = -doY, 

t^o 2 16 

OO 1 1 " 

n > 1 : E(^ + 7;%dk+nY + - ^ dk-A-kY = -d^Y. (4.9) 

fc=0 ^ ^ k=l 



where 

Following our previous notation for the p-spin curves, we use instead the 
parameters 

^™ = E^- (4-11) 

Thus to corresponds to — ti in our notation. There is no to in the expansion 
of Z for the unitary matrix model. Using the Virasoro constraints, we can 
sove for the case in which all the t„, except ti, vanish. The equation (4.9) 
allows one to determine this expansion: 

These terms correspond to the intersection numbers for one, two and three 
marked points. If we set k = —N in our previous result (4.5), in terms 
of Bessel functions, we find that the expression coincides indeed with Y 
given by (4.12). This shows within this expansion the coincidence of the 
two models. Note that the successive terms are independent of A^.The weak 
coupling expansion for general k is given by (4.5). The equation (3.5) was 
missing this k + N, replaced by k ; this is because we were using the replica 
limit to focus on single trace terms. But the true correspondence between 
the unitary and B-model would yield k + N instead of k in (3.5), as one sees 
on (4.5) and we have to set k = —N at the end to recover the unitary model. 



5 Strong coupling expansion 

We now turn to the strong coupling expansion (small Xj) for the B-model. 
We can first use the representation in terms of Bessel functions. Next we 
shall compare with the calculation based on the /-point correlation function 
f/(si, ..., si) for large s and verify that they agree as expected. The ascending 
series for the asymptotic expansion of the modified Bessel function is given 

by 

W-(;^)-E ,„,r(l\L+i) (^■^) 

Putting this expansion into the determinant of (4.2), we find for k = —N for 
the unitary matrix model, 

^ = ^(^ + iV ? ^^ + 2(iV2 - 1) (^ ^^)' ~ 2N{N^ - 1) ^ ^^ 

10 



iV^ — 2 1 

^ 6iV(iV2-l)(iV2-4)^Y^'^ ~ 2(iV2 _ i)(jv2 - 4)^Y^'^ i ^'^ 

where C is a constant term. 

We now consider the strong couphng expansion of the B-model (1.3) 
from the correlation functions f/(si, ..., sj) in its dual representation. Strong 
coupling means small A's. We are thus looking at the large s expansion for 
the Fourier conjugate variable s ~ ^. 

Let us begin with the one point function U{s), given by (3.2). The shift 
M — 7- (m — l)/2 gives 

^^, , 1 / du 4 — /u+ 1\'' 

U(s) ^ -J -e^r^ [—-^) . (5.3) 

As before there is an s-independent constant if we replace the exponential by 

one which does not affect the coefficients of the expansion, and we consider 

thus 

1 ~ 4™ / rfw 1 fu + l\^ . ^, 

^+(^) = 2?;;^/2^R3TfI^ • ^'-"^ 

The contour integral is around the origin, but since the integrand falls off 
at infinity, we can replace the contour by an integral over the discontinuities 
accross the cuts which run from (1, cxd) and (— oo, —1). Then one finds 

(5.5) 

which vanishes as it should when A; is a positive or negative integer. The last 
integral is an Euler beta function: 

We now go to A; = —N as discussed earlier, dropping obviously the van- 
ishing factor sin nk . Note that the Gamma functions give 

which yields only even powers of 1/A^ in a 1/A^-expansion, as it should for 
the initial unitary model. Finally, up to an additional constant and to an 
overall factor we find 

11 



du 1 


fu + l\'' 


2 . , f°^ , (x+lf-"" 

SmTTK / dx- -; 


2m (m2 _ l)m 


U-1/ 


TV Jl (X - 1)^+™ 



This expression agrees, as claimed earlier, with the single trace terms of the 
strong coupling expansion of logZ (5.2) obtained from the determinant of 
Bessel functions. 

The result is an asymtotic expansion, as exhibited by the presence of De 
Wit-'t Hooft poles in A^ [17] for any integer A^, which imply that the 1/A^ 
expansion is certainly only asymptotic. Since the parameter s is a conjugate 
Fourier variable, the use of one single s in U{s), rather than several as in 
U{si ■ ■ ■ si) limits ourselves in the strong coupling (small A) region to terms 
of the form 

j^=j:>^T (5.9) 

If we wanted to compare our results with the full result for the unitary integral 
we would need multiple Si to generate terms of the form J2iX^^ Sj A™^... 
Therefore the expression (4.10) has to be compared with the strong coupling 
expansion of the unitary model in which we would keep only the terms with 
one single trace (Ei A™ = tr{C^Cy"''). 

The strong coupling expansion of the unitary model has been studied long 
ago and we have tables for the expansion due to the work of Samuel [3] who 
studied the model 

ZIj= I dUexp [/3tTAU + /3tTBU^] (5.10) 

where the f3 is the coupling constant. The expansion has the form 

n=l '''• ai,...,«„ 



tiAB trABAB MAB}^ 

N ' ^ N ' ^ N ' 



(5.ii; 



with the restriction ai + 2a2 + ■ ■ ■ + nan = n. (The superscript c in C?^_^^___^^ 
refers to the connected part, the expansion of log Zu rather than Zu). The 
calculations of [3] have been performed up to order /3^°. Without in fact 
losing generality we can substitute (3 A = C, f3B = C^ and use those results 
for the unitary model (1.1). Let us first focus on the coefficients with one 
single trace, namely a„ = 1 and a/ = for / < n. We denote the coefficient 
of Tr(Ai?)" simply by C„ and we find in reference [3] 

1 ^ 

Cl = 1, C2 = —T7^ 7) C'3 



N^-V " (iV2-l)(iV2-4) 

r = ^1 

' (Ar2-l)(iV2-4)(Ar2-9)' 



12 



These numbers agree completely with the coefficients of s^™" in (5.8) and the 
comparison shows that at all orders the result for (5.11) would give for the 
single trace terms. 

If we want to use the same strategy to compare our results with the terms 
of the strong coupling expansion which involve product of several traces (or 
powers of one trace), we have to consider higher point functions. For instance 
the product of two traces requires now to consider the large (si, S2) expansion 
of the two point function f/(si, S2). We know this function explicitely: it is 
given by the double integral (2.9). It has a connected part U{si, S2) which is 
obtained by taking the cycle of length two in the determinant: 

duidu2-, 77 _e"i("i+n) "2(«2+s2) 

[U2 - Ml + S2)[U2 - Ml - Si) 

x(^^^±^)^^^l±^)^ (5.14) 

U2 Ml 

We follow the same strategy : change variables Ui — )■ SiWi, U2 — )■ S2U2 expand 

1 
the exponentials e»"(''+i) in powers of 1/s, but we have also to deal with the 

pole terms ., , \ and Kt-, — 7. There we take a dissymetric limit 

■^ S2(l+M2) — SlUl S2U2 — S\(l+Ui) •' 

in which S2 >■ S\ and replace 

(5.15) 



i 


-'- \ ^ Ml /SlNm 
S2^0(l+«2)'"+l^S2^ 

1 - (X^u^)^ S, 

- Z^ m+1 \ ) 
■52 „=o "2 ^2 


82(1 + M2) - S\U\ 

1 


S2M2-Si(l+Mi) 



(5.16) 

The result is a sum of factorized contour integrals in mi and U2 similar to 
the one that we encountered in (5.4). These integrals may be deformed into 
integrals of discontinuities over the real axis, and after division by the factors 
sin vrA; we can continue to /c = — A^. The integrals of the discontinuities are 
of the type 

N-h 

dx- -^^ = B(N -b-l,b-a-l), (5.17) 

(1+x)^-" ^ ' '' ^ ' 

with the Euler beta function B{p, q) = T{p)T{q)/T{p + q). The coefficient of 
l/s™^l/s™~^ is then given as an explicit sum : 

^ B{N -mi-p2,Pi+P2 + 2mi + 1)B{N -m2+pi + l, 2m2 - Pi - P2 - ^) 

p7^2 (mi +pi +P2 + 1)!("^2 -Pi -P2 - 1)! 

(5.18) 

13 



In the reference [3] we read the first coefficients Ci^m of (^^^)(-jv— ) namely 

Ci,i = 1/{N' - 1) 

Ci_2 =-l2/{N^-l){N^ -A) 

Ci^s =120/{N^-1){N^ -A){N^ -9) 

Ci4 =-lQ80/{N^ -1){N^ -4){N^-9){N^-1Q). (5.19) 

We now compare with the coefficient of the term — ij^rr in the strong couphng 

expansion of f/(si,S2) (5.18). For instance the coefficient of l/siS2 appears 
as a sum of three terms 

N\ 1 1 



2 14 

+ Ar2(Ar2_i)2] = s^siim -l){m - 4)' ^^'^^^ 

note that the sum if a function of iV^ as it should since the unitary integral 
(1.1) has an expansion in even powers of 1/A^. Let us go further and quote 

from [3] the first two 6*2,™. coefficients of { ]\^ ){ ]\^ ) 
„ _ 18{3N^ - 7) 



r =_48n (3^' - 13) .,211 

'•^ {m-iy{m-4){N^-9){N^-i6)' ^' ' 

The sums (5.18) lead to a coefficient of I/sis'^ equal to 
iV2 , 5 



slsl^m{m - l)(iV2 -A){N + 1){N + 3) 
5 



N^m - 1){N^ - 4){N - 1){N - 3) N^{N^-iy{N^-4) 
1 6(3iV2 - 7) 



sfsi(Ar2-l)2(iV2_4)(iV2_9) 

and for l/siS2 

N^ . 42 



(5.22) 



slsYN'^{N'^ - 9)(iV2 _ 4)(jY2 _ i)(Ar + V){N + 4) 

42 

+ iV2(iV2 - 9)(Ar2 _ 4)(jY2 _ Y){N - \){N - 4) 

60 
+ 



Ar2(iv2-i)2(Ar2_4)(Ar2_9) 
1 48(3A^2 _ 13) 

sisf (Ar2 - l)2(Ar2 - 4)(iV2 - 9)(iV2 - 16) ' 
14 



(5.23) 



Again the various sterms of the sums combine and give, as expected, functions 
of N"^ . The comparison with Samuel's results (5.21) is remarkable. Still 
there are combinatorial factors which make the agreement incomplete : for 
instance a factor 10 is missing in the expansion of f/(si,S2) for C2,3 which 
corresponds to the choice 5!/2!3! of assigning si and S2 to a particular trace, 
i.e. the number of ways of splitting 5 objects into two groups of of 2 and 3 ; 
a factor 3 has to multiply our coefficient of 1/ s\s\, i.e. the number of ways of 
splitting the 4 powers in 2 and 2. More generally a factor (mi +m2)!/mi!m2! 
if nil 7^ ^^2 and (2mi)!/2(mi!)^ if irti = 1712- 

Clearly this strategy may be extended to /-point functions (> 3), for any 
/. For instance, if / = 3, we take the connected part of the 3 by 3 determinant 
in (2.9), 

1 



(m2 - Ml + S2){U3 -U2 + S3){ui - M3 + Si) 

+7 T7 T7 V (5-24) 

(M3 - Ml + S3)(M2 -U3 + S2){Ui - M2 + Si) 

and repeat exactly the same expansions with some ordering of Si, 82,83 as 
above, to obtain finally the coefficient of s]~™^s^™^sj™^. For instance, we 
obtain at the lowest order the coefficient of sr^S2"^S3'^ as 

^ ^ (5.25) 



Si82Ss{m-l){N^-4) 



which, up to a sign, coincides with Ci^i^i, the coefficient of {Tt{AB))^. 

We have thus verified that the functions f/(si, ..., s;) provide the gener- 
ating function for the coefficients C^{N) for finite A^. The strong coupling 
expansion, in the presence of the external source C^C has also been analysed 
by Virasoro type recursion relations [5]. However our explicit formulae for 
U{si, ...,8i), from which one may obtain both the strong coupling and the 
weak coupling expansions, are very powerful. 

6 Summary 

We have found that the unitary matrix model may be replaced by a general- 
ized j9-spin higher Airy model which yields the same weak and strong coupling 
expansions. After a duality which replaces this alternative {p, r) model by 
a Gaussian model in an external matrix source, we can compute explicitely 
the Fourier transform of the /-point correlation functions f/(si, ..., s;) (2.9). 
This provides when p = —2 and r = — 1 a generating function for the weak 
and the strong coupling expansions of the unitary matrix model Zu in (1.1). 
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The number / correponds to the total number of traces involving the external 
source matrices CC'^ in the expansions. 

In the weak coupling region, the replica method provides the terms of 
the expansion ordered by the number of independent traces involved in the 
expansion. The zero-replica limit for the one-point function t/(s), determines 
the single trace terms and we have verified the agreement with the weak 
coupling expansion of the unitary matrix model. 

In the strong coupling region, the /-point functions f/(si, ..., S/) remark- 
ably produce the correct results of the strong coupling expansion of the uni- 
tary matrix model with an external source, with its involved iV-dependence, 
De Wit-'t Hooft poles, and allows one to go much beyond the terms which 
had been tabulated long ago. 

Therefore this other matrix model provides finite integrals, on one vari- 
able for the single trace terms, two variables for product of two traces, etc., 
which yield both the weak and the strong coupling regions. It also shows an 
unexpected connexion between the unitary matrix model and the generalized 
Kontsevich model which we know is a generating function for intersection 
numbers [9, 18, 19, 20]. 

The Gross- Witten phase transition between the two regions may also be 
recovered on the integrals for U{s) ; it comes from the logarithmic potential 
of the matrix model (1.3). The applications of these techniques to other 
values of p than —2 with a logarithmic potential may be interesting, and it 
leads also to integral representations for f/(si, ..., s/) as in (2.9). 

Finaly let us note the phase transition in the unitary matrix model 
has been discussed over the last years by several authors in the context of 
string/black hole physics [21, 22, 23]. The reformulation presented here may 
be of interest in this other context as well. 
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